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Abstract 

Wc show how to compute the exact partition function for lattice statistical-mechanical models 
whose Boltzmann weights obey a special "crossing" symmetry. The crossing symmetry equates 
partition functions on different trivalent graphs, allowing a transformation to a graph where the 
partition function is easily computed. The simplest example is counting the number of nets without 
ends on the honeycomb lattice, including a weight per branching. Other examples include an Ising 
model on the Kagome lattice with three-spin interactions, dimers on any graph of corner-sharing 
triangles, and non-crossing loops on the honeycomb lattice, where multiple loops on each edge are 
allowed. We give several methods for obtaining models with this crossing symmetry, one utilizing 
discrete groups and another anyon fusion rules. We also present results indicating that for models 
which deviate slightly from having crossing symmetry, a real-space decimation (renormalization- 
group-like) procedure restores the crossing symmetry. 

1 Introduction 

Over the decades, many techniques have been developed to "solve" classical lattice models., i.e. com- 
pute their partition function in the thermodynamic limit. These typically apply when the Boltzmann 
weights obey the Yang-Baxter equation [Ij. It has become clear that only models with special prop- 
erties are amenable to such treatments; these models are usually called "integrable" . Generic lattice 
models are certainly not integrable, but many of the simplest and most interesting models are, and 
much has been learned about physics as a consequence of these exact results. 

In this paper we develop a very simple technique to compute exact partition functions of some 
classical lattice models. While only a very special class of models is amenable to this technique, such 
models do not seem to have been analyzed before; in general they do not appear to be special cases 
of known integrable models. The degrees of freedom of the models we study live on the edges of any 
lattice (or actually, graph) with trivalent vertices, such as the honeycomb lattice, or equivalently, the 
sites of the Kagome lattice. In several of the simplest and most interesting cases, the models have 
nice geometric interpretations. 

To give an example, the simplest case corresponds to counting the number of "nets" without ends 
on the honeycomb lattice. In this case there are two possible states on each edge, corresponding to 
covering the edge with a net, or leaving it empty. The counting problem becomes non-trivial by the 
requirement that the net has no ends, i.e. no configurations are allowed where only one of the three 
edges touching a vertex is covered. We do allow all three edges touching a vertex to be covered, which 
corresponds to a "branching". A typical net configuration is illustrated in figure [T] below. We show 
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that the number of such net configurations on any planar graph with an even number N of trivalent 
vertices is exactly 

This result can be extended to allow for an arbitrary weighting per branching, as we will explain in 
section [2l 



+ (^1 . (1) 




Figure 1: A typical configuration of closed nets. 

Such computations are possible because in this special class of models, the Boltzmann weights 
obey a very special property we name "crossing" symmetry, by analogy with Feynman diagrams in 
field theory. Crossing symmetry allows us to show the partition functions on different lattices related 
by an "F-move" are identical. Repeating F-moves allows us to transform the lattice into one where 
the partition function is easily computed by diagonalizing an M x M matrix, where M is the number 
of states allowed on each edge of the graph. This procedure of changing the lattice is reminiscent of 
techniques involving the Yang-Baxter equation (in particular that used to solve the eight- vertex model 
on the Kagome lattice [T]), but simpler and much less general. It is also reminiscent of numerical 
techniques used in some classical [2| [3| and quantum two-dimensional models , but our technique 
is not only exact, but can be easily implemented analytically. 

The outline of this paper is as follows: In section [21 we define crossing symmetry and F-moves for 
unoriented edges and vertices. We show how to solve models with this crossing symmetry, illustrating 
this by nets without ends. In section[3]we show how certain loop gases have crossing symmetry, and we 
study some of the properties of these loop models. In sectionS]we solve more general "vertex" models, 
where both the edges and the vertices of the graph can be oriented. One such example generalizes 
the counting of dimer configurations on a Kagome lattice. In section [5] we discuss several methods 
to generate models with crossing symmetry. In section [5T] edges are labeled with the elements of a 
finite group, whereas in section [521 they are labeled with the elements of a quantum theory of anyons 
(i.e., the objects of a modular tensor category, or the primary fields of a rational conformal field 
theory) . In section 15.31 we briefly discuss the relationship of these models based on anyon theories 
to the quantum string-net models of Levin and Wen [5j . Finally in section [6] we consider a real- 
space decimation transformation which is exact for models with crossing symmetry and we give some 
examples of these flows in section 16. li In section 16.21 we consider models which are perturbed away 
from crossing symmetry and we show that under this decimation transformation they "flow" back to 
having crossing symmetry. Finally, in section [7] we summarize our results and discuss further possible 
directions. Several appendices are included to give more details of certain calculations. 
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2 Crossing symmetry and F-moves 



The statistical-mechanical models we study are defined on any lattice or graph Q with trivalent 
vertices, such as the honeycomb lattice. For simplicity we restrict to those Q that can be embedded 
on a sphere, but most of our results will apply to any graph of trivalent vertices. 

The degrees of freedom live on the edges of this lattice, and there are M possible states labeled 
by an index i € . . . M — 1. To allow for more general "vertex" -type models, one can orient each 
edge, but we defer this (only slightly more complicated) case to section [4] below. Here the Boltzmann 
weights depend only on the three states around each vertex, so we label them by w{j, k, I). 




= w{j,k,l) = w{k,lj) = w{k,j,l). 



For now we will further assume that w{j,k,l) = w{k,j,l), but in section H] below we will study the 
more general case where w depends on the order of its arguments. The Boltzmann weight of each 
configuration is the product of all these vertex Boltzmann weights, and the partition function is then 
the sum over all possible configurations. For a model on Q 

Z{G)= Yl n ^Uv,k.,k) (2) 

edge labels vertices v 

where jy, ky, ly are the labels on the edges touching vertex u in a given configuration. Note that we 
have defined the model so that the partition function depends only on the topology of Q and the 
vertex weights w. 




a ^ J 5 

Figure 2: The four types of configurations for a model with M = 2 types of edges. For closed nets, 
there are no endpoints, so /? = 0, while nets have weight per unit length 1 (a = 7 = 1) and branches 
have a weight z (5 = z) . 

A simple example to keep in mind is that of closed nets. Here M = 2, with the two possible 
states on each edge corresponding to the net covering a edge or not. The four types of configurations 
around each vertex and their weights are displayed in figure [2j We use the word "closed" to mean 
that net ends are forbidden, i.e. the Boltzmann weight /3 = 0. There are 8 possible configurations 
at each vertex, the four in the figure with the other four found by rotating the two middle pictures. 
We will show that we can count the number of net configurations without ends, including a weight 
z per branching. The corresponding Boltzmann weights are thus chosen to forbid ends (/3 = 0), to 
give a weight per unit length of 7 = a = 1, and 6 = z. The partition function with these Boltzmann 
weights then reduces to a sum over closed nets C: 

Znets = ) (3) 

C 
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where B is the number of branchings. A typical closed-net configuration is shown in figure [TJ Note 
that, unlike the domain walls in an Ising or Potts model, barbell- type configurations are allowed (i.e., 
the covered edges need not separate two regions from each other). 

The models we solve are invariant under what we call a crossing symmetry. This is a relation 
between a product of the Boltzmann weights of two neighboring vertices 

M-l M-1 

^ w{i,j, m)w{m, k,l) = ^w{l,i, n)w{n, j, k) (4) 

m=0 ra=0 

as illustrated by Fig. [3l 




Figure 3: The definition of crossing symmetry. 



For the model to be solvable by our method, this relation must hold for any choice of i,j,k,l. It is 
simple to check that the Boltzmann weights displayed in figure [2] satisfy crossing symmetry for the 
case of closed nets with weighted branchings. For this case, the only time the sum over the internal 
labels involves more than one term is when all four external legs are covered by nets. 

In general, a useful way of expressing crossing symmetry comes from rewriting the weights w{i,j, k) 
as a set of M symmetric matrices <I>^*) with elements ^^j^ = w{i,j,k). The crossing condition Eq. H] 
then can be rewritten as 

m n 

where here and henceforth all sums run from to M — 1. Crossing symmetry is therefore equivalent 
to the statement that all of the matrices ^'^^^ commute. 

The "crossing" name comes from analogy with field theory, where each trivalent vertex represents 
the fusion of two operators into a third under the operator product expansion. Thus multiplication 
by the matrix <I>*^*) can be thought of as fusing with a particle of type i, such that the matrix element 
means fusing a particle of type i with with one of type j to get a particle of type k. Crossing 
symmetry then means the four-point function can be decomposed into a sum of intermediate states 
using these operator products, with the result independent of the order in which the decomposition 
is done. While our notion of crossing symmetry is much less intricate than that in field theory, we 
do adopt another convenient piece of language from conformal or topological field theory, where such 
decompositions are quite useful: we call the process of exchanging the pictures on the left and right 
of Fig. [3] an F-move. 

Crossing symmetry in our context means that the partition function on graphs related by F-moves 
is identical. Namely when u and v are adjacent vertices, and Fuv{Q) is the graph obtained from Q by 
doing an F-move on u and then 

z{Fug)) = z{g) . 

For this to be true, it is essential that all states be summed over: the equality is obviously not true 
configuration by configuration. 
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The equality of partition functions under crossing symmetry allows us to relate any planar graph 
of trivalent vertices to one where the partition function can be simply calculated. For example, 
consider the following series of F-moves: 






We applied the F-move once in the first step and thrice in the second step, reducing the graph to 
a one-dimensional chain of bubbles. Any planar graph of trivalent vertices can be similarly reduced 
to a chain of bubbles, because each individual face on the graph can be reduced by F-moves to be 
of bubble form. Repeating this for a graph with no dangling edges yields a one-dimensional closed 
chain of bubbles. Since under an F-move the number of vertices in the graph remains the same, given 
N vertices in the initial graph, using the crossing rule we can reduce it to a bubble chain containing 
N/2 bubbles. 

In this form, the model can then be solved trivially by transfer matrix. We write the weight for 
a single bubble, summed over internal lines, as 

n 

3 




m 



The resulting transfer matrix with elements Tjj is an M x M matrix, in terms of which the partition 
function is 

Z{g) = Tr[T^/2] . (6) 

Since the partition function is invariant under F-moves, this applies to any graph comprised of N 
trivalent vertices. Thus computing the partition function of models with crossing symmetry requires 
only diagonalizing an M x M matrix! 

Rewriting the weights in terms of the real symmetric matrices as above, the crossing condition 
means that the matrices all commute. As a consequence, the ^^^^ have the same eigenvectors and can 
be simultaneously diagonalized. Thus we define real orthonormal eigenvectors v^"'^ and eigenvalues 

0(»,a) ^J^g^^ 

k 

In this language the transfer matrix can be written as Tij = Tr[$^*)$(^)] = or as 



m,n 



a,n 



where we used X^a'^l"^'"^^ ~ ^I'm- The eigenvalues of T are thus given by 



a)i2 



(7) 



(8) 



so that the partition function is given by 



a \ n J 



N/2 



(9) 
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For closed nets as shown in Fig. [2] and partition function ([3]) , the transfer matrix is 

T..s=[l 2+%.) , (10) 

which has eigenvalues 

The partition function is thus 

^nets = (A+)^/2 + (A_)^/2 (11) 

The number of such closed net configurations is given by setting z = 1, yielding the result ([U quoted 
in the introduction. The free energy per vertex in the limit of a large number of vertices is ln(A+)/2, 
which for z = 1 is ln[2 cos(7r/10)]. This partition function allows the average number of branchings 
to be computed easily as 

{B) = Z— In(Zncts) • 

az 

A more general two-state model comes by relaxing the restriction that there be no net ends. Such 
a model can also be naturally described as an Ising-type model on any graph Ai of corner- sharing 
triangles. A4 is the medial graph of Q: each vertex of Ai corresponds to an edge of G, and such 
vertices are connected by an edge of Ad if the corresponding edges of Q meet at a vertex of Q. For 
G comprised of trivalent vertices, Ai is comprised of corner-sharing triangles, so for example when 
G the honeycomb lattice, Ai is the Kagome. The degrees of freedom live on the sites of A^, so that 
the configurations in Fig. [5] are rewritten as in Fig. [U where a = w{ ), /3 = w^-l ), etc. 



- + + + 

- - - - - + + + 

a /3 7 (5 

Figure 4: The same weightings of the two state model as shown in Fig. [2] but written in the language 
of an Ising model. 

The crossing relation is not satisfied for arbitrary Boltzmann weights in this two-state model, 
but rather requires that the weights obey. 

a7 + /3<5 = + 72 (12) 

The special case of the closed net model (a = 7 = 1, /? = 0, and 6 = z) indeed satisfies this crossing 
relation. The transfer matrix is then 

/ a2 + 2/32 + ^2 „^ + 2/37 + 75\ 
^ " I a/3 + 2/37 + 7(5 f3^ + 2^^ + 6^ J ^^"^^ 



The partition function is again expressed in the form (jlip . but the explicit form of the eigenvalues 
does not seem to be particularly illuminating here. 

The solution (jl2p of the crossing relations is the most general solution with iVf = 2 states on each 
unoriented edge. We find it remarkable that of the three-parameter family of distinct Boltzmann 
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weights (one weight can always be rescaled to one), a two-parameter family satisfies the crossing 
relation. One can work out the analogous solution for M = 3 states on each edge; we give it explicitly 
in the appendix lA.li There are 10 distinct configurations around each vertex (up to rotations), and 
a three-parameter family of solutions to the crossing relations. Obviously such general solutions will 
get even more unwieldy as the number of allowed states on each edge is increased further. However, 
various special cases can be quite elegant; in the next section, we discuss one of them. 

One thing worth noting is that there is no way to make the two eigenvalues equal in the M = 2 
models with positive Boltzmann weights. (The eigenvalues are equal only in the following three cases: 
/3 = 7 = and Q = (5;a = 7 = and /3 = (5;/3 = (5 = and a = 7.) There is no physically interesting 
phase transition as a function of the Boltzmann weights, a fact that remains true for all models 
solvable via crossing symmetry. This is perhaps not surprising given that the partition function is 
the same as that of a one-dimensional system. A phase transition would occur when the two largest 
eigenvalues of the transfer matrix cross, as this gives a cusp in the free energy. If all the Boltzmann 
weights are non-negative, then all of the entries in the transfer matrix are non-negative, and by the 
Perron-Frobenius theorem we cannot have a degenerate largest eigenvalue. If some Boltzmann weights 
are zero, the only way in which such a degeneracy can occur is if the transfer matrix decomposes into 
blocks. In the one-dimensional model, this type of (first-order) phase transition is somewhat trivial: 
In one phase all of the incoming edge variables are from one set and in the other phase all of the 
incoming edge variables are from a completely disjoint set. (For example, for the above mentioned 
case of /3 = 7 = 0, the transition is between a system that is all spin up to a system that is all spin 
down). However, in more complicated cases with more possible states on each edge, the nature of 
this phase transition is much less obvious. 

3 Loop models 

We showed in the previous section [2] that whenever the Boltzmann weights on a trivalent graph 
satisfy the crossing relation ([4]), the partition function can be easily computed. Our next step is to 
understand better what sorts of models can be solved in this fashion. In this section we describe 
a general family of models having crossing symmetry that are naturally written in terms of loops. 
In the next section we will introduce vertex models and present a several more general classes of 
solutions to the crossing relations. 

In this section we describe models of closed non- intersecting loops on Q. Each of the M states 
i = . . . M — 1 on each edge is pictured by drawing i loops on each edge. The configurations are 
then required to be comprised of loops that do not cross at each vertex. For M = 2, this corresponds 
to setting z = in the closed net model, as apparent in figure [2l Counting these M = 2 loop 
configurations is trivial: the partition function is then simply Z = 2^/^+^, as can be seen directly by 
noting that the loops are domain walls for an Ising model on the dual graph. For M = 3, however, 
the loop model is not trivial: the five types of vertices where loops do not cross are illustrated in 
figure m 

It follows from appendix lA.ll or an easy direct calculation that the weights of our M = 3 loop 
model satisfy the crossing relation for a = 5 = e = 1 and ^ = /i — There thus is a one-parameter 
family of solutions where the partition function can be computed exactly by reducing the graph to 
bubbles by F-moves. In this case the transfer matrix simplifies to 
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Figure 5: The five types of configurations and their weights in the three-state loop model. The middle 
three may be rotated, giving 11 possible configurations at each vertex. 

which has eigenvalues 2(/i^ + 1), 2(/i^ + 1) and 1 + We will see below that this degeneracy of the 
largest eigenvalues, which indicates the system is tuned to a phase transition (for any value of /i), is 
generic for loop models. 

To systematically generalize this solution to higher M, it is useful to recall another context in 
which the crossing relation (j4]) arises. This is as a consistency constraint for the fusion rules for 
operators in a theory of anyons (or equivalently, a rational conformal field theory or topological field 
theory). In such theories the vertex weight w{j,k,l) is non-zero when one of the particles is a "bound 
state" of the other two particles. The crossing constraint then guarantees that a "singlet" bound 
state of four particles can be found by considering any grouping pair by pair. 

A simple case to illustrate this idea is if M — )■ oo and each state j labels the representation of 
the group or algebra SU{2) with spin j/2. Then we let w{j,k,l) = 1 if the tensor product of the 
representations of spin j/2 and spin k/2 contains the representation of spin 1/2, and w{j,k,l) = 
otherwise. Thus w{j, k,l) = 1 when I G \ j — k\/2, \ j — k\/2 + 1, . . . {j + k) /2; this is commonly known 
as the triangle inequality for addition of angular momentum. Note that this definition is symmetric 
in j, k and /. These weights do satisfy the crossing constraint which has a simple interpretation 
in the SU{2) context: If one takes the tensor product of four representations of spin i/2, j/2, k/2 and 
1/2 respectively, either side of the crossing equation corresponds to the number of ways this tensor 
product gives the trivial representation. For example, if all four representations are spin 1/2, each 
side is 2: 



(1/2) (1/2) » (1/2) (1/2) = {(0) + (1)} » {(0) + (1)} = (0) + (1) + (1) + (0) + (1) + (2) . 



The crossing relation shows that to obtain the correct number, it does not matter how the repre- 
sentations are initially paired to do the tensor products. The idea of using a fusion algebra to give 
solutions of the crossing equation turns out to be quite general. In fact, in section [5^2] we will explain 
how the fusion rules of any consistent anyonic theory can be utilized to find a solution of the crossing 
constraint. 

Fusion provides a nice way of finding more general loop models solvable by our method. As above, 
group the vertex weights into a set of M matrices with entries <I>^*^ = w{i,j,k). Satisfying the 

crossing constraints is equivalent to demanding that the matrices ^^^^ all commute with each other. 
For the three-state loop model satisfying the crossing relation, we have ^^'^^ the identity matrix, while 




(14) 



These matrices satisfy 




^$(2) + $(0) 



$(1)$(2) 
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Ignoring the factor of fi, the first of these obviously resembles the fusion rules for SU{2), with fusion 
of the state j akin to taking the tensor product with the spin-j/2 representation (i.e., resembling two 
spin 1/2's fusing to form either spin or spin 1). The second also resembles the SU{2) fusion rules, 
with the caveat that ^^^^ vanishes for this case with M = 3. 

This suggests defining loop models with higher values of M via analogous recursion relations. Let 
= q,q/ with I the M X M identity matrix, and 



$(1) 





"0 
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(15) 



i.e. <I>[.^^ = aiSk,i+i+akSk,i-i with k,l = . . . M—1. This is the most general set of Boltzmann weights 
for fusing one line with k lines to get either k + 1 or k — 1 lines, consistent with non-intersecting loops. 
In the M = 3 case above, = 1 and ai = /i. Then define the remaining weights for fusing more 
than one line via the the matrix recursion relation 



(16) 



The Boltzmann weights defined by the recursion relation ()16p satisfy the crossing constraint 
because all of the matrices <I>('^^ commute with each other. This is easy to prove. Suppose 
commutes with all ^^•'^ for j < i. Using the recursion relation it is easily shown that <I>(*"'"^) also 
commutes with all ^I'^-'^ for j < i + 1. Since ^^^^ is proportional to the identity, it commutes with all 
so by induction all commute. Moreover, w{j, k, I) = as defined here is symmetric in all 
three indices; we prove this in the appendix IB. 2[ 

The w{j, k, I) defined by this recursion relation therefore make sense as Boltzmann weights obeying 
the crossing condition. These loop models therefore can be solved by our method for any M. In 
appendix IB. 31 we show that for even M the eigenvalues of T are always two- fold degenerate, whereas 
for odd M eigenvalues come in pairs except for one remaining unpaired eigenvalue (which cannot 
be the largest eigenvalue of T). As mentioned above, by the Perron- Frobenius theorem, such a 
degeneracy of eigenvalues only occurs if the transfer matrix is block diagonal. Indeed, it is easy to 
see that this is precisely what happens — the transfer matrix preserves the even or odd parity of the 
number of loops on an edge, and is hence block diagonal. 






Figure 6: Illustrating how the recursion relation (|16p results in a non-intersecting loop model. When 
four lines enter from the left, and two enter from the right, the triangle inequality and conservation 
of parity (Eqs. fTTI and llSp then require that the number of lines going out the top must be 2,4 or 6. 

This solvable model defined by (|16p can be represented as a non-intersecting loop model with the 
index A: = 0, . . . , M — 1 of the edge representing k lines running along each edge. This correspondence 
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is illustrated in Fig. [UJ Here, when indices j and k come into two edges of the vertex, a non-vanishing 
weight occurs only if the third vertex has index / with 

\j-k\<l<j + k (17) 
mod (j + /c + /,2) = . (18) 

as required by the above described triangle inequality. Since these two relations hold for k = 1 from 
(I15p . induction using (I16p shows that they remain true for all k. By drawing pictures analogous 
to those in Fig. [6l it is easy to see that these two relations are the requirements for defining a 
non-intersecting loop model, with the M — 1 loops the maximum allowed on each edge. 

If one sets all the aj = 1, all the Boltzmann weights are either or 1. These fusion rules are those 
of the quantum-group algebra Uq{sl2), or the conformal field theory SU{2)m-i- The "truncation" of 
allowing only M loops on each edge is familiar in these contexts. We will explain in more detail in 
section [5^2] how such cases are related to anyon fusion rules. 

4 Vertex models 

In this section we generalize our procedure to solve more general "vertex" models. We still have M 
possible states on each edge of any planar graph with trivalent vertices. Now, however, we orient 
each edge by placing an arrow on it. In vertex models, there is a conjugation operation given by a 
permutation P of [1, ... , M] of order two (i.e., = 1). We denote the conjugate of an index i as i , 
and the conjugation is implemented by reversing the orientation of the edge: 

i _ i 

If all states are self-conjugate (i.e. P is the identity), then we can omit all arrows and the model 
reduces to the unoriented case treated above. The Boltzmann weights of these lattice models depend 
on the three states around each vertex, here taken in counterclockwise order. We define a weight 
function w{i,j,k) = w{j,k,i) = w{k,i,j) corresponding to a trivalent vertex where incoming edges 
have labels and k. 




The partition function is still defined by ([2]) as before, now being careful that in the weights w{i,j, k) 
the edge variables must always be listed in counterclockwise order. 

When we can define an orientation for each vertex (for example if the graph is embedded into an 
orientable surface), then it is possible to have w{i,j, k) ^ w{j, i, k). This is also possible for unoriented 
edges. For M = 2,3 however, crossing symmetry requires w{i,j,k) = 'w{j,i,k); see appendix lAl In 
the remainder of the current section we consider only cases where w{i,j, k) = w{j, i, k), but in section 
15.11 we show how models based on finite groups with larger M do allow weights to depend on the 
orientation of the vertex. 

Crossing symmetry for these oriented vertex models is similar to the unoriented case. Pictorially, 
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m 



n 



'J k 
i.e. in the most general case, the crossing constraint becomes 

y^w{rn,i, j)w{k,l,m) ='y^'w{l,i,n)w{j,k,n) . 



(19) 



The procedure of reducing each graph to a chain of bubble graphs via F-moves is the same as in the 
unoriented case, and the transfer matrix becomes 



Ylr, 



As in the unoriented case, the partition function Z{Q) = Tr[T''^/^] can be rewritten in terms of the 
eigenvalues of the fusion matrix. 

In the simplest vertex model, M = 2 and the two states are conjugate to each other. This is the 
only M = 2 vertex model; if the states are self-conjugate the model reduces to the unoriented case 
considered in section [2j Configurations here can be labeled by putting an arrow on each edge, so that 
the weight of each vertex depends on how many arrows are pointing in and out. 

tt;(in, in, in) = a, w (out, out, out) = 5, 

?i;(in, in, out) = /3, w(out, out, in) = /3 . (20) 

Since every edge appears in two vertex weights, once as in and the other as out, these weights can be 
rescaled as _ _ 

a-f\^a, a-^ X~^a, /3 ^ X/3, /3 ^ (21) 

for any A without changing the partition function. The nontrivial crossing constraint is 

aa = j3j3 . 

The transfer matrix is 

^ _ ( 3/313 + aa 2af3^+ 2/3^ \ 



which using the crossing condition has eigenvalues 4/3/3 it 2( y a/3^ + a(3^). 

A simple but non-trivial special case is /3 = a = 0. In this case there are four non-trivial vertex 
configurations — only all out, or two ins and one out (in three rotations) are allowed. This case corre- 
sponds to counting close-packed dimers on the medial graph Ai , comprised of corner-sharing triangles 
as described in section [2l The vertex configuration in our model with all arrows out corresponds to 
having no dimers on the corresponding triangle, while a vertex with only arrow one out corresponds 
to having a dimer on the edge of the medial graph connecting two in-pointing edges of the original 
graph [6j. Setting 5 = 1 then means that /3 is the weight per dimer. Using our formalism then gives 
_l_2^3/2 £qj. ^YiQ eigenvalues of the transfer matrix, so that the partition function is 



2(2/33/2)^/2 N/2 even 
N/2 odd 



(22) 
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This was derived for Kagome directly from a transfer matrix approach long ago [7]. Here N is the 
number of vertices of the honeycomb graph or triangles of the Kagome. 

The M = 3 vertex model has one state its own conjugate (0 = 0) with the other two conjugates of 
each other (1 = 2). Thus each configuration can be pictured by arrows for 1 and 1, and empty links 
for 0. In appendix El we show that crossing symmetry does not allow for the weights here to depend 
on the orientation of the vertex. This allows for 10 distinct Boltzmann weights at each vertex (up to 
rotations), labeled as in the unoriented case. The general solution for the M = 3 vertex model can 
be found in the appendix \A.2\ again giving six conditions. 

A variety of special examples of the crossing-symmetric M = 3 vertex model can be found by 
demanding that some weights be zero, or that they obey some additional symmetry. An example of 
the former is the analog of the "no net end" condition /3 = 7 = 0, for which one subsequent solution 
is (in the notation of appendix [A]) 

a = 2(/), 6 = (/)Ve, V = 0^'^/e^ M = '/'^C/e^ • 

An example of the latter comes from demanding that the weights be invariant under "charge- 
conjugation" symmetry, i.e. 

w{i,j,k) = w{i,j,k) (23) 

so that here 13 = 'j, 6 = e, rj = and k = fi. The crossing constraint then results in two possible 
solutions for the remaining weights in terms of a, (3 and 6: 

4> = 6, K = C = ^(/32 + 2j2-aJ) , 

or 

(j) = a-S + /3^/6, K-i = P{l-(t)/5), K + i = 5{5 + (t))/P . 
5 More general solutions 

The purpose of this section is to find additional examples of statistical models solvable via cross- 
ing symmetry. We will pursue two routes, one based on using finite groups, and the other using 
anyon theories/modular tensor categories. We also comment on a connection between these models 
and solvable quantum two-dimensional lattice models in section 15.31 but do not attempt a complete 
classification of models solvable by this method. 

5.1 Models based on finite groups 

Here we construct a vertex model for any finite group G, where each directed edge is labeled with 
group elements j, k, . . . G G. Reversing an arrow on an edge corresponds with taking the inverse 
of a group element, so that the conjugate i = . We will show that this model obeys crossing 
symmetry when the Boltzmann weights for a vertex are given by 

i(;(i, j. A;) = l^[conjugacy class(ij/c)] . (24) 

In other words, if the elements entering the vertex (in counterclockwise order) are k, one multiplies 
these elements together to form the group element g = ijk. Two group elements g, h are in the same 
conjugacy class if there exists some element q in the group such that g = qhq~^. Each conjugacy 
class can then be assigned an arbitrary non-negative Boltzmann weight, which we have here denoted 
by W. Since ijk is in the same conjugacy class as kij we must have w{i,j, k) = w{j, k, i) = w{k^ 



12 



as required. We note in passing that models defined this way bear a resemblance to the recent work 
of Ref. [8]. 

For an abelian group, each element is its own conjugacy class, so we can take an independent 
weight for each group element. It is easy to prove that this satisfies the crossing symmetry (jl9p . 
which here amounts to 

^W{ijm)W{m-^kl) = ^W{iln)W{n-^jk) . 

m n 

Letting gi = ijm, g2 = m^^kl, hi = iln, and /i2 = n~^jk, we have gig2 = /11/12 = x. Then the 
crossing relation is equivalent to 

W{gi)W{xigi)-') = ^ Wihi)W{xihi)-'). 

91 hi 

These are obviously the same for any x. 

For a nonabelian group, ijk in general is not equal to jik, so w{i,j,k) 7^ w{j,i,k). For this 
to make sense we must require all the vertices as well as the edges of the graph Q to be given an 
orientation. The proof of crossing symmetry here is not difficult but requires a bit of elementary 
group theory [9]. The character xnid) of a group element g is simply the trace of the matrix D^^\g) 
corresponding to g in the representation R. The character depends only on the conjugacy class of g, 
and any function of a conjugacy class can be expanded in terms of characters xr{9) of irreducible 
representations R. We can thus write 

w{i,j,k) = J2 ^RXniijk) 
R 

where the Aji are arbitrary coefficients, and the sum is over all irreducible representations. The 
required crossing condition, Eq. \T9\ can then be rewritten as 

^ J2 ^R^R' XRiklm)xR'{nr^ij) = X] X] ^R^R' XR{jkn)xR'{n~^li) (25) 

m R,R' n R,R' 

Note that Xr{(^^) = XR.{^o). The following lemma 

V XR{ab)xR' {h-^c) = \^5RR.XR{ac) (26) 

proves that the two sides of (|25p must be equal to each other, thus showing that these models always 
satisfy crossing symmetry. Here ur is the dimension of representation i?, |G| is the order of the group 
G, and the sum over h is over all elements of the group. 

This key lemma is proved from the following matrix manipulations: 

Y,Di'^{ah)DfP{b-^c) = 2:Z?if(a)I)g)(6X)(ri)D(f)(c) (27) 

b b,x,y 

= ^<Ji?i?'5/3-,l)i?(ac) (29) 
hr 
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where the sum on b is over all group elements, and the sums on x, y are over matrix indices. Here 
we have used the fact that D is a representation so that D{ab) = D{a)D{b) as matrices, and in going 
from (j27p to (j28p we have used the well known orthogonality theorem |9J 



sides of (j27p to obtain the desired lemma (j26p . 
5.2 Models based on anyon fusion rules 

In section [3] we found a general class of solutions to the crossing constraints in the unoriented case, 
giving us the sequence of loop models. We noted that these were (when all aj = 1) the fusion rules for 
a theory of anyons. The mathematical structure behind such fusion rules is called a modular tensor 
category^ and is familiar in physics not only in anyon theories, but also in rational conformal field 
theory (RCFT) and topological field theory. In this language, the net models with z = correspond 
to "Fibonacci" fusion rules, whereas the loop models are generally called SU{2)m-i fusion rules. 
In this section we explain how in general any such fusion rules can be used to find solutions of the 
crossing relation. The resulting solutions do not have any free couplings as in the loop models above; 
all weights are non-negative integers. However, as with our earlier analysis, one can then add and/or 
perturb couplings and find still more general solutions to the crossing relation. 

An anyon is a particle in a quantum theory in two spatial dimensions whose statistics generalize 
bosonic and fermionic ones. When an abelian anyon is rotated by 27r around another, the wave 
function changes by a phase. In the non-abelian case, the wavefunction is expressed not as a complex 
scalar, but as a complex vector, describing a state within a degenerate space. Exchange of particles 
can cause a change of state within this space. 

An anyon theory can generally be described by a set of M fundamental particle types, labeled by 
i = . . . M — 1 with being the trivial particle (in the corresponding RCFT, these are called primary 
fields/operators, with the identity operator). A pair of anyons behaves as an anyon as well, and the 
fusion rules describe how the statistics of each pair is described in terms of the fundamental types. 
Precisely, the fusion rules are described by a non- negative integer -^jj^ = Nf^, which gives how many 
times the anyon k appears in the fusion of particles of type i and j, labeled schematically as 



In the RCFT context, a non- vanishing Nf^ says that the field k appears in the operator product 
expansion of i and j. Studying such fusion algebras is a major endeavor; reviews of the rules can 
be found in [10] in the RCFT context, and in appendix E of pTj in the anyon context. A simpler 
overview can be found in [12]. Various examples and a start of a classification scheme can be found 



Any consistent anyon theory /RCFT/modular tensor category defines a classical lattice model 
with crossing symmetry. Each edge label corresponds to an anyon type, and the Boltzmann weight 
w{i,j,k) = Nijk. The fusion rules are symmetric under permutation of particle types: N^f^ is 
independent of the ordering of the indices. Because fusion is associative, such weights automatically 
obey the crossing constraint 




k 



in [13]. 




(30) 



m 



n 
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By our procedure, this gives a solvable model. 

Anyon theories (modular tensor categories), require a great deal of mathematical structure that 
is unnecessary for defining a classical statistical model with crossing symmetry. For example, one 
must have a consistent rule for what happens when two anyons are moved around each other (i.e., 
braided) in two dimensions. A slightly more general mathematical structure which does not define 
such a braiding is known as a unitary fusion category [H]. The objects in such a category also obey 
fusion rules and have the same associativity condition (|30p . It therefore can be used to define a 
classical statistical model with crossing symmetry, even though in general the objects of a unitary 
fusion category cannot be thought of as being actual particles. Even the unitary fusion categories 
have a great deal of mathematical structure that is not necessary for defining a statistical model with 
crossing symmetry (for example, all of the weights are integers). It is an interesting open problem to 
fully classify the most general algebras that would have the necessary crossing symmetry. 

Several well-known examples of fusion rules of anyon models underlie the statistical models de- 
scribed above. The Fibonacci fusion rule describe two types of self-conjugate anyons, which we label 
/ and X. The "identity" particle / fuses according to its name: 1x1 = 1 and I x X = X , while the 
anyon X fuses as X'^ = I+X. The resulting non-zero weights come from A'^m = Nixx = Xxxx = 
with other Nijk not related by permutations vanishing. This gives the net model with z = 1; the 
identity particle corresponds to leaving the edge empty, while the X particle corresponds to covering 
it. So-called Ising fusion rules describe three self-conjugate particles. The particles are the identity 
/, and two additional particles X and Y (often known as a and ip) fusing as X'^ = I + Y, XY = X 
and Y'^ = I with Nyyy vanishing. This is precisely the loop model with M = 3 and ao = ai = 1, 
identifying X and Y with edges having one and two loops respectively. Note that with ai = /i = 1, 
the weight for a vertex with three double lines meeting indeed vanishes, as apparent from (I14p . For 
the loop models with general M and all aj = 1, all A^jjfc are zero or 1, and these correspond to 
SU(2)m-i fusion rules. (Note that Ising and SU{2)2 modular tensor categories are slightly different, 
but they have identical fusion rules.) 

Yet another three-state unoriented (self-conjugate) example is given by the fusion rules [LSI 

X^ = I + Y, XY = X + Y, Y'^ = I + X + Y . (31) 

This case, with Nxxi = Xyyi = XxxY = Xxyy = Xyyy = 1 and others besides permutations 
vanishing, corresponds to having nets with labels X and Y on each strand such that XXX vertices 
are forbidden. As with the loops, we can allow for more general weights. Using the notation of the 
appendix |A] for the Boltzmann weights, the corresponding matrices of weights are 

/a 0\ /O 5 0\ /O e\ 

$(^) =050, =[6 11 f^] , '^^^^ = /i K 

\0 e/ \0 1^ kJ \e k ^/ 

where we have kept f3 = j = (f) = 0to resemble the fusion rules (|3ip . These matrices commute for 
more general couplings than the fusion rules, namely 

a = 5 = e = 1, fi^ + — fj,^ — rjK = 1 , (32) 

giving one constraint in addition to having ^^^^ remain the identity matrix. With these conditions, 

further generalizing the loop recursion relation. It seems very likely that any anyon theory can be 
deformed to generate a larger family of weights that satisfy crossing symmetry. 
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5.3 Duality and the relationship to quantum string-net models 

The models we describe have a great deal in common with the Levin- Wen quantum-mechanical 
"string-net" lattice models [5]. These string-net models were cleverly constructed to give lattice 
realizations of certain topological quantum field theories. We will quote some of the results of this 
work and we refer the reader to the original paper [5j for further details; see also [15] for an alternative 
physical viewpoint. 

As with our statistical models, the Levin- Wen models assign a label i G . . . M— 1 to each edge of a 
trivalent graph. We may choose these labels from an anyon theory (or any unitary tensor category) or 
from a finite group; in the latter case, the resulting Levin- Wen model become in essence a lattice gauge 
theory. There is a one-to-one mapping between edge-label configurations in our lattice model and 
orthonormal basis vectors of the Hilbert space in the string-net model. The precise connection between 
the two occurs when our statistical model consists of vertex weights w which are all either zero or 
one, so that a configuration is either "allowed" (meaning it has weight one) or "disallowed" (meaning 
it has weight zero). In the quantum string-net model, the disallowed configurations correspond to 
a so-called vertex excitation of the model. The Hamiltonian of the model can be tuned to give 
such excitations very high energy and thereby such excitations can effectively be eliminated from 
(low-energy) consideration. We have no need to write the actual Hamiltonian, luckily, as it is quite 
complicated. 

Thus we now consider the quantum string-net models in the absence of any vertex excitations. 
The classical partition function then is the dimension of the remaining low-energy Hilbert space, i.e. 
the number of configurations of the edge variables with only allowed vertices. In the string-net model 
this space can be described in terms of the ground state and all possible excitations which do not 
involve the vertex excitations. The non-vertex excitations are all associated with the faces of the 
model. In particular, each face can can be labeled with a quantum number chosen from the same 
set of M labels (the same labels as we had put on the edges), which in some abstract sense can be 
thought of describing a "magnetic flux" through the face. Moreover, each flux acts like an anyon and 
may fuse with others in different ways, akin to spins adding together to form either singlet or triplet 
representations. Counting of the different ways they may fuse together allows one to determine the 
dimension of the low-energy Hilbert space in another way. 

In order to describe the connection between our statistical model and the quantum string-net 
models more clearly, it is useful to first show that these statistical models have an interesting geometric 
duality. Consider a trivalent graph, and for simplicity assume that it is a planar graph embedded on 
a sphere. Now construct a dual graph in the following way: draw an edge through each face normal to 
the sphere. Using a minimal number of trivalent vertices, connect together all of the endpoints outside 
of the sphere and connect together all of the endpoints inside of the sphere. Such a construction is 
shown explicitly in figure [71 Using the Euler characteristic on a sphere (Faces -|- Vertices — Edges = 2), 
and the fact that all vertices are trivalent so (Edges = (3/2)Faces), it is easy to establish that the 
number of vertices in the dual graph is the same as the number of vertices in the original graph. Thus 
the partition function of the statistical model on the dual graph is identical to the partition function 
on the original graph. 

In the string-net model this dual graph has a very natural interpretation. Each edge through a 
face is labeled with one of the M possible labels, thus describing the flux through each face. These 
fluxes can fuse together in all possible ways on both the inside and the outside of the sphere, and all of 
these possible so-called "fusion channels" describe the full low-energy Hilbert space of the string-net 
model without any vertex excitations. This is equivalent to counting the number of allowed edge-label 
configurations on the original graph, or finding the partition function of the corresponding classical 
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Figure 7: The graph Q has seven black faces on the surface of the sphere (six hexagonal, and one 
very large face encompassing the remainder of the sphere). The dual lattice has an edge (red) drawn 
through each face and hence through the sphere surface. All of the endpoints inside and outside 
of the sphere are connected together using the minimal number of trivalent vertices. Note that the 
number of vertices in the lattice (black) and the dual lattice (red) is identical (ten in each case) . 

model. 

6 Decimation Transform 

In this section we consider a renormalization-group-like decimation transformation in the spirit of 
Ref. [4], in order to analyze the effective long-distance physics. We show that models with crossing 
symmetry remain crossing symmetric under this transformation. We then show that crossing sym- 
metry is stable: if the weights are varied slightly away from crossing symmetry, decimation causes the 
weights to "flow" back to those having crossing symmetry. We also show that these flows are toward 
"infinite temperature", where all Boltzmann weights are the same. 

The decimation transformation proceeds in steps as shown in Fig. [HI The transformation from 




Figure 8: The steps of the RG transformation 

the left to the middle proceeds via multiple F-moves (applied to the circled bonds, for example). The 
transformation from the middle to the right requires turning a triangle (three vertices) into a single 
vertex via a transformation as in Fig. [9l Mathematically, this transformation takes the form 
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Figure 9: Reducing a decorated vertex to a vertex in the decimation 
w{a,b,c) '^w{a,i,j)w{b,k,l)w{cj,k) 

which constitutes the sum over the degrees of freedom of the triangle. 
6.1 The persistence of crossing symmetry 

Crossing symmetry persists exactly under decimation. Namely, after the decimation transformation 
is applied to a crossing-symmetric set of weights, the resulting set of weights will also satisfies crossing 
symmetry. This is easy to show by drawing pictures: 

a & a b 

.b ^ \y 





d 

This is nicely illustrated by the two-state models. For the M 
sectional the decimation transformation is 

13 2a/32 + 2aa/3 + 4/3^/3 



c d/ \ c 
2 vertex model described in 



These still satisfy the crossing constraint aa = (3(3 if the original weights do. To simplify the resulting 
weights, we can use the rescaling transformation (|2ip to rescale (3 to unity. Furthermore, with crossing 
symmetry we can always eliminate a using the constraint. These allow us to obtain the simplified 
decimation transformation 



r/(3 -I- ?7)2 
4/33(1 + 3r/) (3 + r/)/r/ 

where we have defined r] = (3'^ /a. The rj flow has only a single nonnegative attractive fixed point at 
r/ = 1. Using the rescaling (I2ip again, it is then easy to establish that this fixed point corresponds 
to a situation where all four types of vertex have equal weight. This is a sensible outcome: if we 
rescale a huge number of vertices into a single vertex, the total partition function is almost entirely 
independent of the labels on the incoming bonds. One can then think of this as a flow to infinite 
temperature, where all Boltzmann weights are identical. 
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For the M = 2 unoriented model shown in Fig. U the decimation transformation is 

+ 3/3^(a + 7) + 7^ 

P^a^/3 + 6/3^ + 13^ + 2a/37 + + 2/37^ 
7 ^ + 5^7 + 25/37 + 2/3^7 + 07^ + 7^ 

This indeed preserves the crossing condition 07 + [36 = /3^ + 7^. Except for a few singular start- 
ing points where certain weights are zero, the system will always flow to the following form under 
decimation 

(a,/3,7,5) = K{l,x,x^,x^) 
with K and x some constants. This form is stable under decimation, with 

K K^{l + x^f 

for one step of the decimation transformation. Comparing to Fig. 21 we see that each "up" spin has 
fugacity x and K is the remaining weight of the vertex. 

6.2 Flows to crossing symmetry 

We show here that crossing symmetry is stable under this decimation transformation. Namely, a set 
of weights that almost satisfies crossing symmetry flows toward crossing symmetry under decimation. 
Our analysis is very much in the spirit of a real-space renormalization-group procedure. 

Let us consider Boltzmann weights w = Wc-\- e5w, where Wc is a set of weights having crossing 
symmetry, 5w is a deviation from these values and e <C 1 is a small parameter. The partition function 
is written as 

Z{G)= ^ JJ [wcUv,kv,lv) + eS'w{jy,ky,Q] (33) 

edge labels vertices v 

Now let us expand this partition function order by order in e. The first two terms are 

Z{G) = X] n Wc{jv,kv,k) 

\edge labels vertices v 

+ e ^ ^ 6w{jy,kv,ly) Wc{jv,kv,lv) \ + ■ ■ ■ 

\edge labels vertex V vertices v V / 

The first term is the partition function for the model with the crossing symmetric weights Wc- The 
order e term is a sum over terms, where in each term a single vertex is given the perturbation 5w and 
all other vertices are given the crossing symmetric weights Wc- 

It is useful to introduce Zc{G; iai , j ^03 1 • • • , iop), which is the partition function given by not sum- 
ming over the labels on the edges oi , . . . , Op; on these edges the states are fixed to be ia^ , , 1 • • • , iap ■ 
We thus can write the full partition function as 

Z = Zc{g) + e'^ ^ 6w{jv,kv,lv)Zc{g;jv,kv,lv) + ■■■ 

V jv,kv,lv 
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where jy, ky, W a-re the values of the edges incident on vertex V . There is a shght comphcation if the 
graph has a "tadpole"; at the tadpole vertex jv,ky,lv are not all independent. Here we have a sum 
over only two edge variables, say jy, ky, the weight would be 6w{jy, ky,ky), and correspondingly we 
would have Zc{G; jy , ky). The next term in the series, of order e^, involves sums over two vertices 
Vi and V2- It will generically have a sum over the six edges connected to them, except when the 
two vertices are adjacent to each other there will again be fewer independent sums. Typically the 
partially summed partition function Zc{G]iai,ia2jia3, ■ ■ ■ ,iap) is independent of geometry, because 
F- moves on edges other than ai, . . . Up leave it invariant. Thus it is independent of the shape of Q 
except in situations such as having three fixed edge labels ii, 12,^3 all intersect at the same vertex; 
they cannot be geometrically rearranged by F-moves. 



Figure 10: The steps of the decimation transformation with one marked vertex. In the left two panels, 
the cells adjacent to this vertex are numbered for clarity. 

Diagrammatically, let us indicate that a particular vertex takes the value 5w rather than Wc by 
marking it with a dot, as shown on the left of Fig. [101 To lowest order in e we need only sum over 
diagrams each with a single marked vertex. Now consider performing the above described decimation 
transformation. Since the marked vertex does not obey crossing symmetry, we cannot perform an 
F-move on it and we obtain the figure in the middle of Fig. [TOl Let us focus in on the region around 
the marked dot. We can still use F-moves around this region to then transform to the right of Fig. [lOl 
Near the marked point, note in particular that one can still use F-moves around this object as shown 
in Fig. n] 




Figure 11: F-moves still work when an edge is decorated with a tadpole having a marked vertex 

In the last step of the decimation, we would like to to reduce the decorated vertices to single points, 
as shown in Fig. [9l However, recall that the vertices attached to the marked point are not decorated 
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Figure 12: Four objects to be renormalized into a single vertex 



at all. Note that using F-moves we can rewrite the decorated vertex as 




Equivalently, the bubble can be placed on either of the other two legs. We recall that the bubble is 
precisely the transfer matrix T defined above. Let us then define the matrix R = T~^, and notate R 
as a box. Consider now the decorated vertex 



This vertex when connected to a vertex decorated with a triangle (like the left of Fig. [9]) removes the 
triangular decoration and leaves us with only the tadpole. In a picture, 



Recall that we would like to add up all diagrams of the form of the right of Fig. [10] with all possible 
bonds marked. This can be done by using a renormalized vertex where the renormalized vertex is of 
the form of FigfT2l 

We have thus found the decimation transformation explicitly. Each new vertex is comprised of 
the four diagrams in Fig. [T21 Building a hexagonal lattice out of these renormalized vertices precisely 
reproduces the partition function (and the free energy!) to order e. What is remarkable is that, 
due to the relation shown in Fig. [TTl the new vertices obey the crossing relation to order e — even 
though the original ones did not! Were we to consider perturbations to order we would require 
two iterations of the decimation to bring the system to crossing symmetry. 

Thus we conclude that models perturbed slightly from crossing symmetry flow to crossing sym- 
metry under decimation, in the fashion of a real-space renormalization-group flow. This is natural 
when coupled with the observation that crossing-symmetric models themselves flow to having all 
Boltzmann weights equal. Thus in all cases the flow is toward infinite temperature. Models with 
crossing symmetry therefore are in a "high-temperature" phase. 
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7 Summary and Outlook 



We have shown here how crossing symmetry may be used to solve certain classical statistical models 
with ease. Several of these give solutions to simple-to-pose combinatorial problems, such as the 
number of closed nets on any trivalent graph. We described methods of finding more solutions to 
the crossing constraints, indicating a deep connection with the fusion algebra of anyon theories. 
We also argued by using a decimation transformation that such models behave as if they are in a 
high-temperature phase. 

Several key tasks are left for further research. We have explored only a few of the models with 
crossing symmetry; presumably there are other families of particularly interesting special cases. It 
also might be amusing to study these models on higher- genus or higher-dimensional surfaces, which 
are non-trivial to understand for oriented vertices. A better understanding of the appropriate math- 
ematical classification of models with crossing symmetry remains an interesting yet difficult task. It 
also would be interesting to explore if useful quantum models other than the string-net models could 
be constructed from these classical models. Finally, it would be interesting to explore whether there 
are quantum theories where Feynman diagrams may be resummed using the same type of crossing 
symmetries exploited here. 

We would like to thank Jesper Jacobsen for useful comments. S.H.S. would like to thank Fiona 
Burnell, Curt von Keyserlingk, and Karin Erdmann for helpful conversations. The work of P.F. is 
supported by the U.S. National Science Foundation under the grant DMR/MPS1006549, while the 
work of S.H.S. is supported by EPSRC Grant Numbers EP/I032487/1 and EP/I031014/1. 

A Three-state models (M = 3) 

For three states per links, there are ten possible weights at each vertex. We label these as 

u;(000) = a, 'u;(001) = /3, w{002) = 7, u;(011) = 5 w(022) = e, 
w{012) = (l), w{lll)=r], w{222)=^, w{122) = k, w{2ll) = fi. 

We can also allow the weight to depend on the orientation of the vertex, so that ti;(012) = (p is not 
necessarily the same as that of w(021) = cp. However, the crossing symmetry of the three-state model 
then requires 

a{(t>-4>) = f3{(^-4>) = 7(0-0) = <^(0-0) = e(0-0) = v{<P-4>) = ^(0-0) = <4>-4>) = K<P-4>) = 

which then requires (j) = 4>. This is true independent of whether we are discussing the unoriented 
(self-conjugate) model described in lA. II or the vertex (oriented) model described in IA.2I 

A.l Unoriented edges 

Consider the case where edges are all unoriented, i.e., they are all self-conjugate and we draw the 
edge labels without arrows. Given the weights described above and the constraint (j) = (p there are 
then six remaining independent crossing constraints 

—4)^ — K? — + 5e + rjK + = Q 
I3e — ne — ^(j) + + 6k — (pfi = 
—75 + fi6 + f5(j) — (j)r] + (pK — e/j. = 
/32 - r//3 + 52 + 02 _ _ ^ 
13^ — K'j — acp + 5(j) + ecp — fi^ = Q 
7^ - ^7 + + (/)^ - ae - /3k = 
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and the transfer matrix is 

a" + 2/32 + 272 + 52 + £2 + 2</.2 
a/3 + 2(5/3 + 270 + 5ri + eK + 24>^i 
a7 + 2e7 + 2/3(/> + + 2</>«; + 5^x 



al3 + 2(5/3 + 27(/) + (5// + eK + 2</>/i 

/32 + 2(52 + 2(^.2 + 7,2 ^ k2 + 2/u2 

/37 + 2(5(/> + 2e(^ + + r/^ + 2k^ 



07 + 2e7 + 2/3(/) + + 2(/)k + (5^ \ 
;S7 + 2(5(/) + 2e(/) + + W + 2k/U 
72 + 2e2 + 2,^2 + ^2 ^ 2^2 + ^2 y 



A. 2 Vertex model 



Here we take two of the three states (let us choose 1,2) to be conjugate to each other, while state is 
its own conjugate. We can draw 1 and 2 with arrows and as an edge without an arrow. Given the 
weights described above and the constraint (j) = (p, the remaining crossing constraints in the vertex 
model are 

— (/>2 + (5e + T/^ — K// = 
/3e - /ie - 7(?!) + (5^ = 
(/)2 — a(/) + /37 + (5e — /3k — 7^ = 
13'^ - l2P-a6 + 264) - 7^/ = 
7(5 — «;(5 — /3(/) + er? = 
72 - K7 - ae + 2e(t) - f^^ 







and the transfer matrix is 
a2 + 2(/.2 + 4/37 + 2(5e 
a/3 + 20/3 + 27(5 + e/i + (5k + 2(/)/u 
07 + 2(/)7 + 2/3e + (5^ + 2(/)k + efi 



a-f + 2(/)-/ + 2/3e + (5^ + 2(/)k + a/3 + 2<pf3 + 27(5 + e?? + (5k + 2(^^ 

2(/)2 + /37 + 2(5e + r?^ + Sk/x /32 + 2/^2 + 4(5(/) + 2r]K 

72 + 2k2 + 4e(/) + 2^/^ 2(/)2 + /37 + 2(5e + r?^ + Sk/x 



B Some facts about the loop models 

B.l Proving the truncation <|)(^^) = 

Here we show that for M-dimensional matrices <I)(*^) = 0. Let us define the matrix 

"i-i " 
yU) = (_i)i Yla, «>(^) 



.i=0 



The recursion relation (Eq. [T6|) for Y is given by 

y(i+l) ^ _$(l)y(i) 



a,- 



(34) 



and we want to show = 0. 

Note that all the Y's and ^>'s have the same set of eigenvectors which we have labeled v^"'\ Let 
us apply Eq. [3l]to one of these eigevectors, to obtain the recursion 

yi^'+'^ = -'^n#-«|-iyi^-') (35) 

where yn^ is the eigenvalue of Y^^^ corresponding to the eigenvector v^"'^ and <pn is the eigenvalue of 
${i) corresponding to this eigenvector. 

Consider now computing the eigenvalues of we can write the characteristic equation 

/-(/) ao ... \ 

ao ai ... 

ai -(f) 02 



= det 



... aM-3 —4> OiM-2 

V ... aM-2 -<P , 



(36) 
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ai 



Define X^") to be the truncation of this secular matrix to n dimensions. Thus we have 

= 
= 

= 

and so forth. Our characteristic equation is 

detX(^^)=0 (37) 

Now consider calculating the determinant of X^'-'^^'' by method of minors. We obtain the recursion 
relation 

detX(^+i) = -(/.detX(j) - a2_idetX(^-i) 

We recognize this as the same recursion relation as Eq. [35j Thus, if we choose a value of (j) which 
satisfies Eq. [37] (i.e., for 4> one of the eigenvalues (;/)„), then we must also have y^n^^ = as well. Thus 
all the eigenvalues of Y^^^^ are zero, and must be zero as well. 

B.2 Proving that the loop weights are symmetric 

Next we check that w(j, k, I) = ^^^^ is symmetric in all three indices, which is not very obvious from 
the definition. Since the matrices commute and ^^^^ and are symmetric, the recursion relation 
ensures that each matrix is symmetric, so that w{i, k, I) = w{j, I, k). The remaining symmetry follows 
by explicitly plugging in the expression for ^^^'^ in the recursion relation, yielding 

^j^k^^ + ^j-i'^V^ = ^k'^^^l,i + a,^,<P^^l,^, (38) 
ajw{j + l,k,l) + aj-iw{j -l,k,l) = akw{j,k + 1,1) + ak-iw{j,k - 1,1) . (39) 

Consider instead a set of matrices ^''-'^^ defined so that they have entries ^^j^^ = w{k,j, I). Then (j39p 
implies that 

Since this holds for all j and k, this is the same recursion relation as (I38p . This means if the matrices 
= $(0) and ^'^^'> = <&(^), this recursion requires that ^^^^ = for all j, and that the weights 
have the desired symmetry. The last things to prove therefore are that w{j,0,k) = w{0,j,k) and 
w{j, 1, k) = w{l,j, k) For the first, assume that w(m, 0, k) = w{0, m, k) = ao5m,k for all m < j. Then 
(I39p with / = and the last two entries exchanged implies that 

w{j + 1,0, k) = — {-aj-i6j^i^k + 0!kSj^k+i + ak-iSj^k-i) = aoSj,k-i = w{0,j + l,k) 



a 



Since w{l,0,k) = ao(5fc,i by definition, recursion therefore implies w{j,0,k) = 'w{0,j,k) for all k as 
desired. The proof that w{j, 1, k) = w{l,j, k) follows analogously by setting / = 1 in (f39l) . 
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B.3 Pairing of Eigenvalues 

Consider the matrix Its eigenvalue equation is 

ak+iVk+i + Ok-iVk-i = (p^^^Vk 

with eigenvector with indices k = 0, ... M — 1 and eigenvalue (p^^^ (having defined here a_i = 
om = 0). From the eigenvector we can construct another eigenvector given by 

ik = {-l)''vk 

which has eigenvalue (f)^^^ = —(f)^^\ Thus eigenvalues of ^^^^ must come in pairs. The exception to 
this is if 0^^^ = which occurs if M is odd which corresponds an eigenvector with Vk = for all odd 
k. 

Correspondingly let us define (p^^^ to be the eigenvalue of ^^^'^ corresponding to the eigenvector 
Vk, and let 4'^^^ to be the eigenvalue of ^^^^ corresponding to the eigenvector Vk- From the recursion 
Eq. [16] we see that these eigenvalues also come in pairs 

^(j) ^ (_i)i0(i) 

Thus using Eq. [8] we see that the eigenvalues of T must also come in pairs, with the exception of 
one lone eigenvalue which occurs only in the case of M odd, which corresponds to the eigenvector 
with zero eigenvalue for In the case of M odd, we should ask whether the lone eigenvalue could 
possibly be the largest eigenvalue of T. In fact, this can never occur by the Perron- Frobenius theorem, 
which guarantees that the eigenvector corresponding to the largest eigenvalue must have all positive 
entries. 
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